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DEFORMATION  OF  A  LIQUID  SURFACE  BY  AN  IMPINGING  GAS  JET* 


J1  AN  MARC  \  AND!  N-HROI  C  K 


Abstract.  1  he  deformation  of  a  liquid  surface  due  to  an  impinging  two-dimensional  jet  is  considered 
assuming  potential  How  This  problem  was  solved  numcncallv  for  small  valuer  *>t  the  Fruude  numkr  a  bv 
Olmstead  md  Ravnor  -  1^64.  Bv  using  a  different  numerical  procedure  we  solve  the  problem  for  larger 
values  of  a  up  to  A,  -6.  We  show  that  the  profile  of  the  liquid  surface  contains  a  train  of  waves  far  awav 
from  the  impinging  jet 


C 


■L. 


1.  Introduction.  In  an  interesting  theoretical  and  numerical  work,  Olmstead  and 
Raynor  ( 1964)  examined  the  impact  of  a  vertical  jet  on  a  stationary  horizontal  liquid 
surface  (see  Fig  1 ).  They  derived  the  integro-differential  system  1 1 .1 H  1 .3)  below  for 
the  complex  velocity  fir )  on  the  free  surface  of  the  liquid.  They  simplified  it  by  replacing 
sin  0  by  ft  in  1 1.2)  and  then  solved  the  resulting  system  numerically.  Their  numerical 
scheme  converges,  but  only  with  few  mesh  points  and  for  relativel>  small  va  es  of  A, 
with  indication  of  a  rapid  growth  in  the  error  as  A  increases. 

We  shall  derive  a  numerical  scheme  different  from  theirs  to  solve  (1  1  i-t  1 . 3 1  isee 
§  2».  It  enables  us  to  compute  the  solution  for  arbitrarv  values  of  A  up  to  A  ■■=  A  -6. 
Our  computations  show  that  the  liquid  surface  contains  a  train  of  waves  far  from  the 
jet  (see  Figs.  1  and  2).  These  results  are  discussed  in  $3 

The  problem  is  characterized  by  the  velocity  \  ,  and  the  width  b  of  the  jet  at  >'  ~  x . 
the  gas  density  pc„  the  liquid  density  pt  and  the  acceleration  of  gravity  g.  Taking  b/  n 
as  the  unit  length  and  \  V  as  the  unit  velocity,  Olmstead  anti  Raynor  <  1 964 1  derived  the 
follow  ing  nonlinear  integro-differential  equation  for  the  omplex  velocitv  c  on  the  free 
surface  CJ:: 


,  >  1  f  I  1  -  S *  i  V  » 
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on  CJ^. 


Here  A  is  the  Froude  number  defined  by 

(1.4)  A  ~2 p^Vj/p,  gb. 


The  variable  r  in  ( 1 .  l)-(  1.3)  is  related  to  the  potential  function  <b  on  CJi  by  the  formula 
(1.5)  e  2*  =  (l~r)/2,  ref -1, +11 


2.  Numerical  procedure.  To  solve  (J.1H  1.3)  Olmstead  and  Raynor  (1964)  intro¬ 
duced  the  N  mesh  points  r]  defined  by 


(2.1) 


Ti  =  1 


/  -  I ,  • ,  N. 
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Fl<».  1  Profile  of  the  jet  and  the  liquid  surface  for  A  -  0.6  The  axes  and  a  sketch  of  the  flo*  are  also  shots  n 


Fu.  Z  Profile  <o  the  liquid  surface  for  \ 


Approximating  sin  ft  by  ft  in  M.2»  and  discretizing  tl.li  and  il.2i.  they  obtained  a 
system  of  S  equations  with  tX  unknowns  77  =  r(r; ),  /  =  1,  •  ,  S.  This  system  was 

solved  by  Newton's  method,  using  the  solution  corresponding  to  A  -  0  as  the  initial 
guess.  Converged  solutions  of  the  algebraic  system  of  equations  for  ,V  =  20  and  A  ^  ] 
were  presented  in  their  paper.  In  addition,  they  gave  an  estimate  of  the  error  based  on 
calculation  of  the  vertical  force  on  the  liquid  surface.  It  showed  that  the  error  grows 
rapidly  with  A,  reaching  8%  for  A  =  1.  A  small  part  of  this  error  can  be  attributed  to 
the  approximation  sin  ft  —  6.  How'ever,  the  main  error  in  their  solutions  comes  from 
the  fact  that  their  liquid  profiles  do  not  contain  waves.  Numerical  experiments  show 
that  Olmstead  and  Raynor's  numerical  procedures  always  give  solutions  of  the  discret¬ 
ized  equations  without  waves.  However,  for  any  A  >0,  there  exists  a  critical  number 
of  mesh  points  Nc{\  )  such  that  the  scheme  converges  for  X  <  S,  (A  )  and  diverges  for 
;V  >  S\  (A  ).  In  addition  the  function  /V,  (A  )  decreases  as  A  increases. 

Equations  (1.1  Hi. 3)  are  reminiscent  of  the  integro-ditferential  equation  derived 
by  Vanden-Broeck  and  l  uck  (1977)  and  Vanden-Brocck,  Schwartz  and  Tuck  (1978) 
to  describe  the  flow  behind  a  moving  body.  Using  the  analogy  between  these  two 
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problems  and  following  the  general  philosophy  of  Vanden-Broeck,  Schwartz  and 
l  uck's  method,  the  following  successful  numerical  scheme  was  derived 

First  we  use  (1.5)  to  rewrite  U.lMl.3)  in  terms  of  the  new  independent  variable 
< b •  Next  we  introduce  the  N  mesh  points  </»j  given  by 

(2.2)  <£/  =  (/  - 1  )F,  7  =  1,  ■*,Ar. 

Here  E  is  the  interval  of  discretization.  We  also  introduce  the  .V  corresponding 
unknowns, 

(2.3)  =  0[r\4>t)l  /  =  ],- 

By  symmetry  6X  =  0,  so  only  N  -  1  of  the  0/  are  unknown.  We  shall  also  use  the  N  -  1 
intermediate  mesh  points  0/  ♦  i  /  2  given  by 

(2.4)  <*;•  1/’ =  £(<*>/ +<*/.■>.  /  =  !,-•  •  ,N  -l. 

We  now  compute 

r/.  1,2  =  r[r(<£, . ,  2l] 

in  terms  of  0/  by  applying  the  trapezoidal  rule  to  the  integral  in  ( 1 . 1 1  rewritten  with  the 
new  variable,  with  the  mesh  points  <£/.  The  symmetry  of  the  discretization  enables  us 
to  compute  the  Cauchy  principal  value  as  if  it  were  an  ordinary  integral  The  error 
inherent  in  approximating  the  integral  by  an  integral  over  a  finite  interval  was  found 
to  be  negligible  for  NE  large  enough.  Then  from  r/.i  2  we  compute  r*(d>/.i  : >  and 
from  0/  we  compute  0/.j  2  In  the  discretization  we  have  used  five-point  difference 
formulas  and  four-point  interpolation  formulas  and  obtain  the  results  in  terms  of 
Next  we  substitute  into  (1.2)  the  expressions  so  obtained  for  7/.  1, r*(d>/ .  1  :)  and 
0/ .  1  2  at  the  iV  2  points  d>j  M /  =■  2,  •  •  • ,  S  -  1 .  Thus  we  obtain  a  system  of  .V  2 
nonlinear  algebraic  equations  involving  the  N  -  1  unknowns  0/,  J  =  2.  •  *  *  ,  K.  T  he  last 
equation  is  obtained  by  expressing  0:  in  terms  of  0*,  04  and  0s  by  a  three-point  Fagrange 
extrapolation  formula. 

For  each  value  of  A,  the  A7  -  1  equations  are  solved  by  Newton's  method  w  ith  0/  -  0 
as  the  initial  guess.  The  remaining  part  of  the  computations  follows  closely  Olmstead 
and  Raynor's  work.  The  scheme  was  found  to  be  rapidly  convergent  and  a  solution  of 
the  algebraic  equations  with  an  error  less  than  10  "  was  obtained  in  a  few  iterations 
For  each  value  of  A  and  ;VE,  the  value  of  E  was  decreased,  and  correspondingly  \  w  as 
increased,  until  the  computed  profiles  remained  unchanged,  at  least  to  graphical 
accuracy.  The  procedure  was  then  repeated  for  large  values  of  A E  until  the  results 
became  independent  of  A IF.  Most  of  the  results  presented  in  the  next  section  were 
obtained  with  E  =  0.3  and  N  -  100. 

3.  Discussion  of  the  results.  Typical  profiles  or  the  liquid  surface  fur  A  -  O.b  and 
A  5  are  presented  in  Figures  1  and  2.  In  both  figures  the  horizontal  scale  has  been 
shrunk  to  show'  clearly  the  train  of  uniform  waves  far  from  the  jet.  It  is  convenient  to 
characterize  these  waves  by  their  steepness  s  defined  as 


Mere  Ay  is  the  difference  of  ordinate  between  a  crest  and  a  trough  and  /  is  the 
wavelength.  The  values  of  s  are  presented  as  a  function  of  A  in  Figure  3.  As  A  tends 
to  zero,  the  steepness  .v  tends  to  zero  and  the  waves  can  be  approximated  by  the  sine 
waves  of  linear  theory  (see  Fig.  1).  It  can  easily  bo  verified  that  the  wavelength  /  of 
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F iCi .  3.  Values  of  the  steepness  s  as  a  function  of  the  Froiuie  numbc  A. 


the  waves  satisfies  the  dispersion  relation 

(3.21  L  -  A7t“  tanh  ~  asA-+0. 


As  A  increases,  the  steepness  5  increases  and  the  waves  start  tv)  develop  sharp 
troughs  and  broad  crests  (see  Fig.  2).  These  w^ves  behave  qualitatively  like  gravity 
surface  waves  of  finite  amplitude  turned  upside-down.  This  could  be  expected  since 
the  complete  problem  can  be  interpreted  as  a  negative-gravity  water-wave  problem 
turned  upside-down  (Tuck  (1975)).  For  steep  waves  the  velocity  is  small  at  the  troughs 
and  large  at  the  crests.  Since  the  mesh  points  correspond  to  equal  increments  in  the 
velocity  potential,  they  are  quite  dense  near  the  crests  and  sparse  in  the  vicinity  of  the 
troughs.  This  nonuniform  spacing  of  the  mesh  points  for  A  large  limits  the  accuracy  of 
the  present  numerical  scheme.  Similar  difficulties  w'ere  encountered  before  by  Schw  arts 
and  Vanden-Broeck  (1979)  and  Chen  and  Saffman  (1980).  Accurate  solutions  for 
A  >At  could  not  be  computed  even  with  N  =  120, 
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